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The deuteron elastic structure function A(Q 2 ) has been extracted in the range 0.7 < Q 2 < 6.0 
(GeV/c) 2 from cross section measurements of elastic electron-deuteron scattering in coincidence 
using the Hall A Facility of Jefferson Laboratory. The data are compared to theoretical models based 
on the impulse approximation with the inclusion of meson-exchange currents, and to predictions of 
quark dimensional scaling and perturbative quantum chromodynamics. 
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Electron scattering from the deuteron has long been a 
crucial tool in understanding the internal structure and 
dynamics of the nuclear two-body system. In particu- 
lar, the deuteron electromagnetic form factors, measured 
in elastic scattering, offer unique opportunities to test 
models of the short-range nucleon-nucleon interaction, 
meson-exchange currents and isobaric configurations as 
well as the possible influence of explicit quark degrees of 
freedom §, §. 

The cross section for elastic electron-deuteron (e-d) 
scattering is described by the Rosenbluth formula: 



da 



A(Q 2 ) + B(Q 2 ) tan 2 (-) 



(1) 



where a M = a 2 E' cos 2 (<9/2)/[4£; 3 sin 4 (6>/2)] is the Mott 
cross section. Here E and E' are the incident and scat- 
tered electron energies, is the electron scattering an- 
gle, Q 2 = AEE' sin 2 (6'/2) is the four-momentum transfer- 
squared and a is the fine structure constant. The elas- 
tic electric and magnetic structure functions A(Q 2 ) and 
B{Q 2 ) are given in terms of the charge, quadrupole and 
magnetic form factors Fc{Q 2 ), Fq(Q 2 ), Fm(Q 2 )- 



A{Q 2 ) = F 2 (Q 2 ) + °-r 2 F 2 (Q 2 ) + ^F 2 U (Q 2 ) (2) 



B(Q 2 ) 



-(l+r)F 2 M (Q 2 ) (3) 



where r = Q 2 /4M i 2 , with Md being the deuteron mass. 

The interaction between the electron and the deuteron 
is mediated by the exchange of a virtual photon. In the 
non-relativistic impulse approximation (IA), where the 
photon is assumed to interact with one of the two nu- 
cleons in the deuteron, the deuteron form factors are 
described in terms of the deuteron wave function and 
the electromagnetic form factors of the nucleons. The- 
oretical calculations based on the IA approach Q| using 
various nucleon-nucleon potentials and parametrizations 
of the nucleon form factors generally underestimate the 
existing A(Q 2 ) data ||, f§, ||, ||- Recent relativistic 
impulse approximation (RIA) calculations improve or 
worsen the agreement with the data depending on their 
particular assumptions. There are two RIA approaches: 
manifestly covariant calculations JtJ, || and light-front 
dynamics §, (To). 

It is well known that the form factors of the deuteron 
are very sensitive to the presence of meson-exchange cur- 
rents (MEC) |lj. There have been numerous extensive 
studies and calculations augmenting both the IA and 
RIA approaches with the inclusion of MEC. The prin- 
cipal uncertainties in these calculations are the poorly 
known value of the pnj coupling constant and the pttj 
vertex form factor. Some calculations show also sensi- 
tivity to possible presence of isobar configurations in the 
deuteron pu. The inclusion of MEC brings the theory 
into better agreement with the existing data. 



It is widely recognized that the underlying quark-gluon 
dynamics cannot be ignored at distances much less than 
the nucleon size. This has led to the formulation of so- 
called hybrid quark models that try to simultane- 
ously incorporate the quark- and gluon-exchange mecha- 
nism at short distances and the meson-exchange mecha- 
nism at long and intermediate distances. When the inter- 
nucleon separation is smaller than ~ 1 fm, the deuteron is 
treated as a six-quark configuration with a certain prob- 
ability that results in an additional contribution to the 
deuteron form factors. 

At sufficiently large momentum transfers the form 
factors are expected to be calculable in terms of only 
quarks and gluons within the framework of quantum 
chromodynamics (QCD). The first attempt at a quark- 
gluon description of the deuteron form factors was 
based on quark-dimensional scaling (QDS) The 
underlying dynamical mechanism during e-d scatter- 
ing is the rescattering of the constituent quarks via 
the exchange of hard gluons. The Q 2 dependence 
of this process is then predicted by simply counting 
the number of gluon propagators (5), which implies 
that y/A(Q 2 ) ~ (Q 2 )~ 5 . This prediction was later 
substantiated in the framework of perturbative QCD 
(pQCD), where it was shown |l4| that in leading-order: 

Jim = K(Q 2 )/g 2 ] 5 E m ,„rf™„[in(g 2 /A 2 )]^-^ 

where a s (Q 2 ) and A are the QCD strong coupling 
constant and scale parameter, and j m ,n an d d mn are 
QCD anomalous dimensions and constants. The exist- 
ing SLAC A(Q 2 ) data Q exhibit some evidence of this 
asymptotic fall-off for Q 2 > 2 (GeV/c) 2 . 

The unique features of the Continuous Electron Beam 
Accelerator and Hall A Facilities of the Jefferson Labora- 
tory (JLab) offered the opportunity to extend the kine- 
matical range of A(Q 2 ) and to resolve inconsistencies in 
previous data sets from different laboratories by mea- 
suring the elastic e-d cross section for 0.7 < Q 2 < 6.0 
(GeV/c) 2 . Electron beams of 100% duty factor were 
scattered off a liquid deuterium target in Hall A. Scat- 
tered electrons were detected in the electron High Resolu- 
tion Spectrometer (HRSE). To suppress backgrounds and 
separate elastic from inelastic processes, recoil deuterons 
were detected in coincidence with the scattered electrons 
in the hadron HRS (HRSH). Elastic electron-proton (e- 
p) scattering in coincidence was used to calibrate this 
double-arm system. A schematic of the Hall A Facility 
as used in this experiment is shown in Figure 1. 

The incident beam energy was varied between 3.2 and 
4.4 GeV. The beam intensity, 5 to 120 /iA, was moni- 
tored using two resonant cavity beam current monitors 
(BCM) upstream of the target system. The two cavities 
were frequently calibrated against a parametric current 
transformer monitor (Unser monitor) |15|. The beam was 
rastered on the target in both horizontal and vertical di- 
rections at high frequency and its position was monitored 
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FIG. 1. Plan view of the Hall A Facility of JLab as used 
in this experiment. Shown are the beam monitoring devices, 
the cryo-target, the two magnetically identical spectrometers 
(consisting of quadrupoles Qi, Q2, Q3, and dipole D), and 
the detector packages. 

Scattered electrons were detected in HRSE used in its 
standard configuration consisting of two planes of plas- 
tic scintillators to form an "electron" trigger, a pair of 
drift chambers for electron track reconstruction, and a 
gas threshold Ccrcnkov counter and a segmented lead- 
glass calorimeter for electron identification. Recoil nuclei 
were detected in HRSH using a subset of its detection sys- 
tem: two planes of scintillators to form a "recoil" trigger 
and a pair of drift chambers for recoil track reconstruc- 
tion. The efficiencies of the calorimeter and Cercnkov 
counter were ~ 99.5%, and of scintillators and tracking 
almost 100% for both spectrometers. Event triggers con- 
sisted of electron-recoil coincidences and of a prescaled 
sample of electron and recoil single- arm triggers. 

Electron events were identified on the basis of a mini- 
mal pulse height in the Cerenkov counter and an energy 
deposited in the calorimeter consistent with the momen- 
tum determined from the drift chamber track. Coinci- 
dence events were identified using the relative time-of- 
ffight (TOF) between the electron and recoil triggers. 



Contributions from the target cell end-caps and random 
coincidences were negligible. Elastic e-p scattering was 
measured for each e-d elastic kinematics. The e-p kine- 
matics was chosen to match the electron-recoil solid an- 
gle Jacobian for the corresponding e-d kinematics. Data 
were taken with and without acceptance defining colli- 
mators in front of the spectrometers. 

The elastic e-p and e-d cross sections were calculated 
using: 

da_ _ N ep{d) C ef f 

dQ NiNt(An) M cF 1 ' 

where N ep ^ is the number of e-p(e-d) elastic events, Ni 
is the number of incident electrons, N t is the number of 
target nuclei/cm 2 , (AO)mC is the effective double-arm 
acceptance from a Monte Carlo simulation, F is the por- 
tion of radiative corrections that depends only on Q 2 and 
T (1.088 and 1.092, on average, for e-p and e-d elastic re- 
spectively) and C e ff = CdetCcdtCmi- Here C de t is the 
electron and recoil detector and trigger inefficiency cor- 
rection (2.6%), C c dt is the computer dead-time correction 
(typically 10% for e-d elastic), and C rn i is the correction 
for losses of recoil nuclei due to nuclear interactions in the 
target (0.7-1.8% for protons and 2.8-5.1% for deuterons). 
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FIG. 2. The deuteron elastic structure function A(Q 2 ) 
from this experiment compared to RIA theoretical calcula- 
tions @|. Also shown are previous SLAC data [Q. 

The effective double-arm acceptance was evaluated 
with a Monte Carlo computer program that simulated 
elastic e-p and e-d scattering under identical conditions 
as our measurements. The program tracked scattered 
electrons and recoil nuclei from the target to the detec- 
tors through the two HRS's using optical models based 
on magnetic measurements of the quadrupole and dipole 
elements and on position surveys of collimation systems, 
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magnets and vacuum apertures. The effects from ion- 
ization energy losses and multiple scattering in the tar- 
get and vacuum windows were taken into account for 
both electrons and recoil nuclei. Bremsstrahlung radia- 
tion losses for both incident and scattered electrons in 
the target and vacuum windows as well as internal ra- 
diative effects were also taken into account. Details on 
this simulation method can be found in Ref . (JT8) . Monte 
Carlo simulated spectra of scattered electrons and recoil 
nuclei were found to be in very good agreement with ex- 
perimentally measured spectra. 
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FIG. 3. The present A(Q 2 ) data compared with overlap- 
ping data from CEA ||, SLAC @LBonn §, Saclay §, and 
IA+MEC theoretical calculations pi . 



The e-p elastic cross sections measured with the ac- 
ceptance defining collimators were found to agree within 
0.3%, on average, with values calculated using a recent 
fit to world data of the proton form factors. The e-p 
clastic cross sections measured without the collimators 
were, on average, 2.6% higher than the ones measured 
with collimators. All e-d cross section data taken with- 
out collimators have been normalized by 2.6%. 

Values for A{Q 2 ) were extracted from the measured 
e-d cross sections under the assumption that B(Q 2 ) does 
not contribute to the cross section (supported by the ex- 
isting B(Q 2 ) data [^6)). The extracted A(Q 2 ) values are 
presented in Fig. 2 together with previous SLAC data 
M and theoretical calculations. The error bars represent 
statistical and systematic uncertainties added in quadra- 
ture. The statistical error ranged from ±1% to ±28%. 
The systematic error has been estimated to be ±5.9% and 
is dominated by the uncertainty in (ACI)mc (±3.6%). 
Each of the two highest Q 2 points represents the average 
of two measurements with different beam energies (4.0 
and 4.4 GeV). Tables of numbers are given in Ref. Jl9| . 
It is apparent that our data agree very well with the 



SLAC data in the range of overlap and exhibit a smooth 
fall-off with Q 2 . 

The double dot-dashed and dot-dashed curves in Fig. 
2 represent the RIA calculations of Van Orden, Devine 
and Gross (VDG) § and Hummel and Tjon (HT) [§, 
respectively. The VDG curve is based on a relativisti- 
cally covariant calculation that uses the Gross equation 
p(j| and assumes that the virtual photon is absorbed by 
an off-mass-shell nucleon or a nucleon that is on-mass- 
shell right before or after the interaction. The HT curve 
is based on a one-boson-exchange quasipotential approx- 
imation of the Bethe-Salpeter equation pi] where the 
two nucleons are treated symmetrically by putting them 
equally off their mass-shell with zero relative energy. In 
both cases the RIA appears to be lower than the data. 
Both groups have augmented their models by including 
the pirj MEC contribution. The magnitude of this con- 
tribution depends on the pirj coupling constant and ver- 
tex form factor choices f22]| . The VDG model (dashed 
curve) uses a pnj form factor from a covariant separa- 
ble quark model j?3) . The HT model (dotted curve) uses 
a Vector Dominance Model. The difference in the two 
models is indicative of the size of theoretical uncertain- 
ties. Although our data favor the VDG calculations, a 
complete test of the RIA±MEC framework will require 
improved and/or extended measurements of the nucleon 
form factors and of the deuteron B(Q 2 ), planned at JLab. 
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FIG. 4. The deuteron form factor F d (Q 2 ) times (Q 2 ) 5 (top) 
and the reduced deuteron form factor fd{Q 2 ) (bottom) from 
this experiment and from SLAC Q. The curve is the asymp- 
totic pQCD prediction of Ref. (II) for A=100 MeV, arbitrarily 
normalized to the data at Q 2 = 4 (GeV/c) 2 . 

Figure 3 shows our data in the "low" Q 2 range where 
they overlap with data from other laboratories. The pre- 
vious measurements tend to show two long-standing di- 
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verging trends, one supported by the SLAC data and the 
other one by the CEA || and Bonn [|| data. Our data 
agree with the Saclay data || and confirm the trend of 
the SLAC data. It should be noted that another JLab 
experiment has measured A(Q 2 ) in the Q 2 range 0.7 to 
1.8 (GeV/c) 2 [ pi| . The two curves are from a recent 
non-relativistic IA calculation |2^] using the Argonne v\s 
potential without (dot-dashed curve) and with (dashed 
curve) MEC, and exhibit clearly the necessity of MEC 
inclusion also in the non-relativistic IA. 

Figure 4 (top) shows values for the "deuteron form 
factor" F d (Q 2 ) = y/ A(Q 2 ) multiplied by (Q 2 ) 5 . It is ev- 
ident that our data exhibit a behavior consistent with 
the power law of QDS and pQCD. Figure 4 (bottom) 
shows values for the "reduced" deuteron form factor |2g] 
f d {Q 2 ) ee F d {Q 2 )/F 2 (Q 2 /4) where the two powers of 
the nucleon form factor F N (Q 2 ) = (1 + Q 2 /0.71)~ 2 re- 
move in a minimal and approximate way the effects of 
nucleon compositeness |2q| . Our f d (Q 2 ) data appear to 
follow, for Q 2 > 2 (GeV/c) 2 , the asymptotic Q 2 predic- 
tion of pQCD 0: f d (Q 2 ) ~ MQ 2 )/Q 2 ][ln(Q 2 /A 2 )]- r . 
Here T = -(2C F /5/3), where C F = (n 2 - l)/2n c , f3 = 
11 — (2/3)rif, with n c = 3 and nf = 2 being the numbers 
of QCD colors and effective flavors. Although several 
authors have questioned the validity of QDS and pQCD 
at the momentum transfers of this experiment p7| , p8| , 
similar scaling behavior has been reported in deuteron 
photodisintegration cross sections at moderate photon 
energies p9| . 

In summary, we have measured the elastic structure 
function A(Q 2 ) of the deuteron up to large momentum 
transfers. The results have clarified inconsistencies in 
previous data sets at low Q 2 . The high luminosity and 
unique capabilities of the JLab facilities enabled mea- 
surements of record low cross sections (the average cross 
section for Q 2 = 6 (GeV/c) 2 is ~ 2 x 1CT 41 cm 2 /sr) that 
allowed extraction of values of A(Q 2 ) lower by one or- 
der of magnitude than achieved at SLAC. The precision 
of our data will provide severe constraints on theoretical 
calculations of the electromagnetic structure of the two- 
body nuclear system. Calculations based on the relativis- 
tic impulse approximation augmented by meson-exhangc 
currents are consistent with the present data. The results 
are also indicative of a scaling behavior consistent with 
predictions of dimensional quark scaling and perturba- 
tive QCD. Future measurements, at higher Q 2 , of A(Q 2 ) 
and B(Q 2 ) as well as of the form factors of the helium 
isotopes would be critical for testing the validity of the 
apparent scaling behavior. 
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